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In view of the imminent start of the LHC experimental programme, we use the available indirect experimental 
and cosmological information to estimate the likely range of parameters of the constrained minimal supersymmet- 
ric extension of the Standard Model (CMSSM), using a Markov-chain Monte Carlo (MCMC) technique to sample 
the parameter space. The 95% confidence-level area in the (mo, mi/ 2 ) plane of the CMSSM lies largely within the 
region that could be explored with 1 fb 1 of integrated luminosity at 14 TeV, and much of the 68% confidence-level 
area lies within the region that could be explored with 50 pb _1 of integrated luminosity at 10 TeV. A same-sign 
dilepton signal could well be visible in most of the 68% confidence-level area with 1 fb -1 of integrated luminosity 
at 14 TeV. We discuss the sensitivities of the preferred ranges to variations in the most relevant indirect experi- 
mental and cosmological constraints and also to deviations from the universality of the supersymmetry-breaking 
contributions to the masses of the Higgs bosons. 
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1. Introduction and topical to make the best possible assessment 

___ of the likelihood that the LHC will indeed dis- 

Supersymmetry (SUSY) PEE] is one of the CQver STJSY; based Qn the begt available experi . 

most highly favoured extensions of the Standard phenom enological and cosmological infor- 

Model (SM), and is often considered to be a prime mation ^ Mogt of the current constraints on pos _ 

candidate for discovery at the Large Hadron Col- sMc pnysicg beyond the gM &re negative; in thc 

hder (LHC). With the start of experiments at senge that they refiect the agreement of data with 

the LHC now becoming imminent, it is natural the gMj and get Qnly lower limits Qn the possible 
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masses of supersymmetric particles [4]. Exam- 
ples are direct constraints such as lower limits on 
specific sparticles, e.g., the chargino, and indirect 
constraints such as the lower limit on the possible 
mass of a SM-like Higgs boson |5l6j . 

However, there are two observational con- 
straints that, within the context of SUSY, may 
be used also to set upper limits on the pos- 
sible masses of supersymmetric particles, since 
they correspond to measurements that cannot 
be explained by the SM alone. These hints for 
new physics are the anomalous magnetic moment 
of the muon, (g — 2) M , which appears to dif- 
fer by over three standard deviations from the 
best SM calculation based on low-energy e + e~ 
data |7I8I9I10I11I12IT3] . and the density of cold 
dark matter, f^cDM P3], which has no possible 
origin within the SM. Each of these discrepancies 
has many possible interpretations, of which SUSY 
is just one. Nevertheless, given the strong moti- 
vations for SUSY, which include the naturalness 
of the mass hierarchy and grand unification, as 
well as the existence of a plausible candidate for 
the astrophysical cold dark matter, it is natural 
to ask what (g — 2) M and f^cDM may imply for 
the parameters of supersymmetric models. Any 
such analysis should also take into account the 
constraints imposed by precision measurements of 
electroweak observables (EWPO) and £?-physics 
observables (BPO) such as BR(6 — > 57), where 
most observables agree quite well with the SM. 

In this paper we revisit the indirect information 
on supersymmetric model parameters obtainable 
in the light of these experimental, phenomcnologi- 
cal and cosmological constraints, using a Markov- 
chain Monte Carlo (MCMC) approach, see, e.g., 
Ref. 15 and references therein. This is practi- 
cal only in simplified versions of the minimal su- 
persymmetric extension of the Standard Model 
(MSSM) , in which some universality relations are 
imposed on the soft SUSY-breaking parameters. 
Initially, we work in the framework of the con- 
strained MSSM (CMSSM), in which the scalar 
and gaugino mass parameters, mo, TO1/2, and 
the trilinear coupling Aq are each assumed to 
be equal at the input GUT scale. Furthermore 
as low-energy parameter we have tan/3, the ra- 
tio of the two vacuum expectation values of the 



two Higgs doublets. At the end we also com- 
ment on the possible changes in our results if 
the common soft SUSY-breaking contribution to 
the Higgs scalar masses-squared, rn 2 H , is allowed 
to differ from those of the squarks and sleptons, 
the single-parameter non-universal Higgs model 
or NUHM1. 

There have been many previous studies of 
the CMSSM parameter space [161171181191 
I20l2ll22l23l24l25l26l27l28l29l30l3ll32l33l34l351 
I36I37I38I39I40I41I42I43I44I45I46I47I48] . includ- 
ing estimates of the sparticle masses, and 
a number of these have used MCMC tech- 
niques [34|35|36|37|38|39|40|41|42|43] . These 
have been used to extract the preferred values for 
the CMSSM parameters using low-energy preci- 
sion data, bounds from astrophysical observables 
and flavour-related observables. These analyses 
differ in the precision observables that have been 
considered, the level of sophistication of the the- 
ory predictions that have been used, and the way 
the statistical analysis has been performed. 

Here we use the MCMC technique to sam- 
ple efficiently the SUSY parameter space, and 
thereby construct the \ 2 probability function, 
P(x 2 , Adof). This accounts correctly for the num- 
ber of degrees of freedom, N^oi, and thus repre- 
sents a quantitative measure for the quality-of-fit. 
Hence -P(x 2 ,Ad f) can be used to estimate the 
absolute probability with which the CMSSM de- 
scribes the experimental data. Our probabilistic 
treatment is explained in detail in Sec. [2] 

Many previous analyses found evidence for 
a relatively low SUSY mass scale in the stau- 
coannihilation region, e.g., [32|47|49] . and a mild 
preference for tan/3 ~ 10 was found in [4T]. A 
comparison of Bayesian analyses yielding vary- 
ing results under different assumptions was made 
in (48) . Some differences between analyses may 
also be traceable to the treatments of the BR(6 — > 
57) and (g — 2) M measurements, which we discuss 
in some detail below. 

Our main objectives in this paper are threefold. 
One is to discuss explicitly the prospects for dis- 
covering sparticles in early LHC running, another 
is to discuss the robustness of the fit results by 
analyzing the implications of relaxing the con- 
straints due to (g — 2)^, BR(fr — > S7), Ocdm and 
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other observables, and the third is to discuss the 
extension of the CMSSM results to the NUHM1, 
in which an extra parameter is introduced that 
allows a common degree of non-universality for 
the two Higgs multiplets. 

We find that the 95% C.L. area in the 
(toi/ 2 ,too) plane of the CMSSM lies largely 
within the region that could be explored with 
1 fhr 1 of integrated LHC luminosity at 14 TeV 
in a single experiment, and that much of the 
68% C.L. area lies within the region that could 
be explored with 50 pb" 1 of integrated lumi- 
nosity at 10 TeV (the projected initial LHC 
collision energy). A same-sign dilepton signal 
could well be visible in the 68% C.L. area with 
1 fb _1 of integrated luminosity at 14 TeV, and 
the lightest Higgs boson might also be detectable 
in squark decays with 2 fb _1 of integrated lu- 
minosity at 14 TeV. We find that removing the 
Ocdm constraint has little effect on the preferred 
regions of the CMSSM parameter space in the 
(mo, TO1/2), (tan/3, TO1/2)) and (tan/3, mo) planes, 
apart from expanding the range of mo, particu- 
larly for tan (3 ~ 10. On the other hand, rescaling 
the present error in — may have quite an im- 
portant effect: the preferred ranges in m 1 / 2 and 
mo would expand quite significantly if the error 
on the present experimental discrepancy with the 
SM were to be increased. Conversely, if this error 
could be reduced, e.g., by a more precise measure- 
ment of (g — 2)n and/or a more refined theoretical 
estimate within the SM, the predictions for spar- 
ticle masses could be significantly improved. We 
also discuss the effects of possible variations in 
the errors in BR(6 —> 57) and other observables. 
Finally, we show that our results would not be 
greatly changed in the NUHM1: we leave a more 
complete study of the NUHM1 and the NUHM2 
(in which the masses of the two Higgs multiplets 
are independently non-universal) for future work. 

2. Multi-parameter Fit to Experimental 
Observables 

Important observables used in our analysis are 
listed in Tab. [TJ Some of the EWPO that are 
included in the analysis have not been listed in 
the Table, because they did not change since the 



analysis carried out in [47]; their details can be 
found there. 

The deviation of (<? — 2) M from the SM pre- 
diction by more than 3 a can be easily accommo- 
dated within the (C)MSSM by choosing appropri- 
ately the sign of the Higgs supermultiplet mixing 
parameter,^: sign(/i) = sign(a° xp — o? M ). Conse- 
quently, we analyze in detail the case /j, > 0, and 
discuss the fi < case only briefly. 

The central value of the BR(6 — > sj) con- 
straint has changed slightly because of new ex- 
perimental results: the data/SM ratio in Tab. Q] 
corresponds to the HFAG average BR(6 — > sj)= 
(3.52±0.24) x 10~ 5 [66] and to the NNLO SM cal- 
culation, BR(6 -> S7)= (3.15 ± 0.23) x 10" 5 \5T\ 
(both values refer to the inclusive rate with E 7 > 
1.6 GeV). Despite some interesting recent at- 
tempts to improve the SM prediction of BR(6 — > 
57) (see, e.g., Refs. |85|86|87(88] and references 
therein), following Ref. [55] we still consider the 
above NNLO value as the most reliable SM esti- 
mate. As compared to [47], we have reduced the 
additional theoretical error in the calculation of 
the SUSY contribution, for the following reasons. 
First, it should be noted that all non-perturbative 
uncertainties cancel out in the SUSY/SM ratio. 
Secondly, data force the deviations from the SM 
to be small in BR(6 -> 37), so the SUSY/SM ra- 
tio can be computed to a relatively high degree 
of accuracy 0- A conservative 15% error on the 
b — > S7 SUSY amplitude corresponds to less than 
5% in the BRSUSY/ br sm^ ratio in the region 

where this does not deviate from unity by more 
than 30%. The BR(6 — > sj) constraint, as well 
as all the other flavour-physics constraints listed 
in Tab. [1] have been implemented using the code 
developed in Refs. [731 74) . This includes the lead- 
ing NLO QCD corrections to the supersymmetric 



1 There are two exceptional cases where the theoretical un- 
certainties of the SUSY amplitude can be large: i) the 
SUSY amplitude is about twice the SM one (SUSY/SM~ 
—2) yielding a BR(6 — -> sf) rate close to the SM value 
ii) the overall SUSY contribution is small because of can- 
cellations among independent large terms. Case i) is ex- 
cluded by the B -» X s l+£- constraints |89I90| that we 
take into account in our numerical analysis. We deal 
with case ii) by implementing in our code the leading 
NLO SUSY contributions, that are known within the MFV 
framework 91 1621 . 
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Observable 


Th. Source 


Ex. Source 


Constraint 


Add. Th. Una 


m w [GeV] 


[5tJI5T] 


m 


80.399 ±0.025 


0.010 




|8l53l54l55] 


|7|l0|56| 


(30.2 ±8.8) x 10~ 10 


2.0 x 10~ 10 


m h [GeV] 


[57 58 59 60^ 


[5.6j 


> 114.4 (see text) 


3.0 


RR cxp /RR SM 


Ul UZ UO U^t UOJ 




i.ii ( rn u.u ( Doxp =c u.uozth(SM) 


u.uou 


^\ 7~\ 

m t [^evj 


|oUlol] 


[DYj 


1 70 /I _1_ 1 O 






[68169170] 


dU 


0.1099 ±0.0062 


0.012 


BR(£? S — > 


[71172173174] 


[66] 


< 4.7 x 10~ 8 


0.02 x 10~ 8 


BR^ T JBR™ TV 


73 74 75J 


76 77 78J 


1.15±0.40 [cxp+th] 




BR(i?d — * 


[71 72 73 74 


[66] 


< 2.3 x 10~ 8 


0.01 x 10~ 9 


RT? CXp /T3TJ SM 


m 


[66180] 


0.99 ±0.32 




T^TjCxp mnSM 


[73175] 


ED 


1.008 ±0.014 [cxp+th] 




RT? CXp /T3TJ SM 


[82] 


[83] 


< 4.5 




AM^ p /AAfSM 


[82] 


m 


l.ll±0.01oxp±0.32 th(SM ) 




(AM° P /AM™) 
(AMJ7/AM|M) 


[71172173174] 


|66|84] 


1.09±0.01 eX p±0.16 th(SM ) 




Ae£ p /Ae| M 


[82] 


[84] 


0.92 ± 0.14 [cxp+th] 





Table 1 



List of experimental constraints used in this work in addition to the electroweak observables listed in [47j . 
The top part of the table shows observables that are very sensitive to the MSSM parameter space, the 
middle part lists observables with updated measurements compared to [37] while the bottom part lists 
additional experimental constraints. The values and errors shown are the current best understanding of 
these constraints. The rightmost column displays additional theoretical uncertainties taken into account 
when implementing these constraints in the MSSM. 



contributions [62j and a complete resummation 
of all the relevant large tan/3 effects beyond the 
lowest order 63 64 65J. More recent public codes 
for the evaluation of BR(6 -s- sj) in the CMSSM 
have been presented in Ref. [92193] , 

A significant _B-physics constraint arises also 
from BR(_B — > tv), which represents a powerful 
probe of the (m#±,tan/3) plane [73|74j . However, 
at present both experimental and theoretical un- 
certainties prevent us from fully exploiting the 
potential sensitivity of this observable. In par- 
ticular, the SM prediction suffers from the uncer- 
tainties in the determination of the CKM element 
\V u b\ and of the decay constant fs- Concern- 
ing \V u b\, we use the current HFAG average [66] 
(from combined exclusive and inclusive semilep- 



tonic B decays), while for fs we use the lattice 
result of 78J. An alternative way to reduce the 
theoretical error associated to BR(£> — > tv)sm 
is to consider the ratio BK(B — > tv)/ AMs d , 
where fs drops out and \V u b\ is replaced by 
\V ub /V t d\ = sin 13/ sin 7 [73|74j . However, as 
the experimental error on BR(£> — > tv) is the 
dominant uncertainty, this alternative way does 
not lead to a significant reduction in the er- 
ror. Moreover, it complicates the analysis since 
BR(B — » tv) and AMs d are affected by in- 
dependent SUSY contributions. For these rea- 
sons, we treat the two constraints separately 
More precisely, we treat separately BR(£> — > tv), 
AM Bs , and the ratio AM Bi /AM Bi . Similarly to 
BR(& — > 57), all flavour-physics constraints apart 
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from BR(_B S —¥ are implemented nor- 

malising the observables to the corresponding SM 
values. 

The direct experimental limit on the Higgs- 
boson mass in the SM obtained at LEP [5] is 
m h > 114.4 GeV at the 95% C.L. The corre- 
sponding bound within the MSSM could in prin- 
ciple be substantially lower, due to a reduced ZZh 
coupling or due to different, more complicated 
decay modes of the Higgs bosons [5J. However, 
it has been shown [94195) that these mechanisms 
cannot be realised within the CMSSM, and hence 
the experimental lower bound of 114.4 GeV can 
be applied 0. For our fit we use the full likelihood 
information of the exclusion bound, given by the 
CL s (rrih) value, which is convoluted with a theory 
error on the evaluation of mh of 3 GeV [57], ac- 
cording to the detailed prescription found in [29] . 

The numerical evaluation has been performed 
with the MasterCode that consistently combines 
the codes responsible for RGE running, for which 
we use SoftSUSY 96] , and the various low-energy 
observables. At the electroweak scale we have 
included the following codes: FeynHiggs [57|58[ 
159160] for the evaluation of the Higgs masses and 
a^ USY ; a code based on [731 74j and Superlso [53] 
for the flavour observables; a code based on [50] 
I5T] for the electroweak precision observables; 
MicrOMEGAs [68|69|70] and DarkSUSY [97|98] for 
the observables related to dark matter. We 
made extensive use of the SUSY Les Houches Ac- 
cord [99] in the combination of the various codes 
within the MasterCode. 

The CMSSM parameter space has been sam- 
pled using the MCMC technique. We treat mi/2, 
mo, Aq and tan/3 as free parameters, and the 
Higgs mixing parameter /x and the pseudoscalar 
Higgs mass toa as dependent parameters deter- 
mined by the electroweak vacuum conditions. 

A global x 2 function is defined, which combines 



2 Following Ref. [32] , for simplicity we use this bound also 
in our NUHM1 analysis. As discussed below, the best-fit 
NUHM1 point we find yields mj, well above this bound. 



all calculations with experimental constraints: 

N in V>\i /fobs ffit \2 

2 _ \ ~> (yi ~ fjj \ ~» WSMj JSMj) 

X -2f (T (c i )* + ( r(P i )* + 2r' ^(/smJ 2 

(1) 

Here N is the number of observables stud- 
ied, C'i represents an experimentally measured 
value (constraint) and each Pi defines a CMSSM 
parameter-dependent prediction for the corre- 
sponding constraint. The three SM parameters 
fsM = {Aahad, TOt, mz} are included as fit pa- 
rameters and constrained to be within their cur- 
rent experimental resolution <t(/sm)- 

As indicated in Section [TJ the sensitivity of the 
global fit to different constraint scenarios is stud- 
ied below by removing the SIcdm constraint or 
rescaling the (<? — 2)„ and other experimental un- 
certainties. Since each new scenario represents a 
new x 2 function which must be minimized, multi- 
ple re-samplings of the full multi-dimensional pa- 
rameter space are, in principle, required to deter- 
mine the most probable fit regions for each sce- 
nario and would be computationally too expen- 
sive. 

To avoid this difficulty, we analyze the effect of 
removing the Ocdm constraint by exploiting the 
fact that independent % 2 functions are additive 
and result in a well-defined x 2 probability. Hence, 
a "loose" x 2 function, X 2 00 se> ^ s defined in which 
the term representing the Ocdm constraint is re- 
moved from the original x 2 ■ The X 2 oose function 
represents the likelihood that a particular set of 
model parameter values is compatible with a sub- 
set of the experimental data constraints, without 
any experimental knowledge of Ocdm ■ 

An exhaustive, and computationally expensive, 
25 million point pre-sampling of the Xioose f unc ~ 
tion in the full multi-dimensional model param- 
eter space is then performed using an MCMC. 
The result of this pre-sampling identifies fit re- 
gions which are generally excluded by the consid- 
ered sub-set of experimental data. Any regions 
excluded by the less constrained fit will also be 
excluded with the inclusion of additional exper- 
imental constraints and, in particular, with dif- 
ferent scenarios for the Ocdm constraint. Hence, 
without loss of generality, this pre-sampling pro- 
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cedure reduces the hyper-volume of parameter 
space which needs to be searched multiple times 
over in the context of different constraint scenar- 
ios to a computationally manageable level. 

Constraint terms representing the different 
!]cdm scenarios are then re-instated to form dif- 
ferent x 2 = xfoosc + Xsconario functions, one for 
each scenario studied. The precise values of the 
most probable fit parameters are determined via a 
full MINUIT pOT] minimization of the X 2 for each 
different scenario, but are performed only within 
the general parameter space regions not already 
excluded from the pre-sampling of the Xioose fhnc- 
tion. An MCMC final sampling is subsequently 
used to determine the 68% and 95% likelihood 
contours for each scenario constraint studied. 

Additionally, later on we vary the uncertain- 
ties of ricDM and other constraints using similar 
techniques. This allows us to study and compare 
the effects of such variations on the \ 2 fit and the 
most probable parameters. 

3. Results 

3.1. Predictions for LHC discoveries 

In Fig. [T] we display the best-fit value and the 
68% and 95% likelihood contours for the CMSSM 
(mo, mxn) plane, obtained as described in Sect. [2] 
from a fit taking into account all experimental 
constraints listed in Tab. [1] as well as the con- 
straints from the additional electroweak observ- 
ables listed in [37]. We also show in the up- 
per panel of Fig. [T] various LHC sparticle dis- 
covery contours for 1 fb -1 of good-quality data 
in a single experiment at a centre-of-mass en- 
ergy of 14 TeV. The ATLAS and CMS collab- 
orations have each published 5-cr discovery con- 
tours in the CMSSM (mya, mo) plane for Ao = 
and tan/3 = 10 |101|102|103"] . Their contours are 
generally very similar, and the solid brown con- 
tour displayed is that published by CMS for the 
most sensitive jets + missing Et search. This 
contour is insensitive to Aq, which affects primar- 
ily the third-generation sparticle masses, since the 
main discovery channels involve gluinos and first- 
generation squarks. The discovery contours are 
also not very sensitive to tan/3, since the gluino 
mass is insensitive to this variable, and the first- 



generation squark masses are also not very sen- 
sitive to tan p. Therefore, it is a reasonable first 
approximation to compare our 68% and 95% like- 
lihood contours directly with the discovery con- 
tours given for A$ = and tan p = 10 fixed, 
particularly since the best fit has a similar value 
of tan (3. 

The parameters of the best-fit CMSSM point 
are JB1/2 = 310 GeV, m = 60 GeV, A = 
240 GeV, tan/3 = 11 and /1 = 380 GeV 0, 
yielding the overall xV N dof = 20.4/19 (37.3% 
probability) and = 113.2 GeV 0. The over- 
all value of the \ 2 a ^ the minimum is some- 
what pushed up by the value of m^, which is 
uncomfortably low. However, it is acceptable 
within the higher-order calculational uncertain- 
ties expected in the FeynHiggs code that we use 
here, Sm h hc ° « 3 GeV [57] . As we discuss be- 
low, this slight tension is removed in the NUHM1 
model, which correspondingly has a somewhat 
lower overall \ 2 (yielding a similar fit probabil- 
ity for the two models). The spectrum at the 
best-fit CMSSM point is shown in the left panel 
of Fig. [2] It is interesting to note that the best-fit 
CMSSM point and the corresponding spectrum 
are quite similar to the well-known SPSla bench- 
mark point [104] , whose phenomenology at future 
colliders has been studied in considerable detail 
(see, e.g., |105|106|107] ). 

Comparing the 95% likelihood contour pro- 
vided by the multi-parameter fit with the 1 fb _1 
LHC discovery contour, we see that the former 
is almost entirely contained within the latter, im- 
plying that, if the CMSSM were correct, the LHC 
would be almost 'guaranteed', with 95% confi- 
dence, to discover SUSY with 1 fb _1 of good- 
quality data at 14 TeV. We also display in the 
upper panel of Fig. [1] contours representing the 
5 a discovery reach with 1 fb _1 at 14 TeV for 
4-jet events with and without a charged lepton, 
for same-sign dileptons |102j . denoted by SS, and 
(with 2 fb" 1 ) for the lightest MSSM Higgs boson 



3 Here and later, we quote CMSSM and NUHM1 input 
mass parameters with 10 GeV accuracy. 
4 The CMSSM fit quality has improved relative to [47] pri- 
marily because of the new value of rat and the inclusion 
of more observables, that are generally highly consistent 
with the CMSSM. 
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Figure 1. The (mcm^) plane in the CMSSM for tan/3 = 10 and Aq = 0. The dark shaded area at low 
mo and high TO1/2 is excluded due to a scalar tau LSP, the light shaded areas at low wii/2 do not exhibit 
electroweak symmetry breaking. The nearly horizontal line at m 1 / 2 ~ 160 GeV in the lower panel has 
m~± = 103 GeV, and the area below is excluded by LEP searches. Just above this contour at low mo 
in the lower panel is the region that is excluded by trilepton searches at the Tevatron. Shown in both 
plots are the best-fit point, indicated by a filled circle, and the 68 (95)% C.L. contours from our fit as 
dark grey/blue (light grey/red) overlays, scanned over all tan/? and Aq values. Upper plot: Some 5tr 
discovery contours at ATLAS and CMS with 1 fb _1 at 14 TeV, and the contour for the 5cr discovery 
of the Higgs boson in sparticle decays with 2 fb _1 at 14 TeV in CMS. Lower plot: The 5 a discovery 
contours for jet + missing E T events at CMS with 1 fb" 1 at 14 TeV, 100 pb" 1 at 14 TeV and 50 pb" 1 
at 10 TeV centre-of-mass energy. 
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Figure 2. The spectra at the best-fit points: left — in the CMSSM with m = 60 GeV, m 1/2 = 310 GeV, 
A = 240 GeV, tan/3 = 11, and right — in the NUHM1 with m = 100 GeV, m 1/2 = 240 GeV, 
A = -930 GeV, tan/3 = 7, m 2 H = -6.9 x 10 5 GeV 2 and /x = 870 GeV. 



produced in cascade decays of sparticlcs [102 , de- 
noted by h. We see that the same-sign dilepton 
discovery region largely covers the 68% likelihood 
region of the CMSSM (mo,m 1 / 2 ) plane. Thus, 
this signature could serve as a clean signal capa- 
ble of confirming the supcrsymmetric interpreta- 
tion of any jets + missing Et signal observed in 
initial LHC running. On the other hand, the re- 
gion where the lightest Higgs boson could be dis- 
covered in cascade decays of squarks with 2 fb _1 
at 14 TeV lies largely between the 95% and 68% 
C.L. contours. 

We have used PR0SPIN02 [108] to estimate the 
variation of the discovery reach of the LHC jets 
+ missing Et search as a function of the inte- 
grated luminosity and the centre-of-mass energy. 
We display in the lower panel of Fig. [1] (green) 
dot-dashed and (red) dashed contours represent- 
ing, respectively, the discovery reaches expected 
with 100 pb" 1 at 14 TeV and 50 pb" 1 at 10 TeV. 
We see that the 68% likelihood contour is well 
covered by the 14 TeV/100 pb _1 discovery reach, 
and even the 10 TeV/50 pb _1 reach would be suf- 
ficient to discover SUSY at the best-fit point, in- 



dicated by a filled circle in Fig.[T] The lower panel 
of Fig. [1] also displays the regions of the CMSSM 
(mojirii^) plane that are excluded by chargino 
searches at LEP and by sparticle searches at the 
Tevatron [4|109|110j . The region excluded by the 
LEP Higgs search is sensitive to tan (3 and > 
is subject to theoretical uncertainties, and, more- 
over, the experimental Higgs likelihood function 
is not a simple step function. Hence, it is not 
shown in Fig. [T][B 

3.2. Sensitivity to experimental con- 
straints 

The above analysis assumed the default imple- 
mentations of the experimental, phenomenologi- 
cal and cosmological constraints discussed in the 
previous Section. We now discuss the possible ef- 
fects of relaxing (or strengthening) some of the 
key constraints, starting with the relic cold dark 
matter density, SIcdm- 

It is well-known that this constraint essentially 
reduces the dimensionality of the MSSM param- 

5 However, for orientation, we note that if mo = 0, tan f3 = 
10 and Aq = the evaluation with FeynHiggs yields a 
nominal value of ra/, = 114.4 GeV for rrii/2 =307 GeV. 
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eter space by one unit, fixing one combination of 
the parameters with an accuracy of a few %. For 
example, in the CMSSM for any pair of fixed val- 
ues of Aq and tan /3, the Ocdm constraint largely 
determines too as a function of TO1/2, except for a 
discrete ambiguity associated with the coannihi- 
lation strip, the focus-point strip and the rapid- 
annihilation funnel that appears at large tan/?. 
Therefore, one might expect that dropping the 
Ocdm constraint would have a strong effect on 
the preferred region of the CMSSM (m 1 / 2 , mo) 
plane shown in Fig. [T] 

There are various possible reasons why one 
might consider dropping the dark matter con- 
straint. Perhaps the neutralino is not the LSP? 
Perhaps i?-parity is not quite conserved? Perhaps 
the early thermal history of the Universe differed 
from that usually assumed when calculating the 
relic LSP density? Perhaps Nature is described 
by some generalization of the CMSSM such as a 
model with non-universal SUSY-breaking contri- 
butions to the Higgs scalar masses (NUHM), in 
which case values of too very different from those 
in the CMSSM might be permitted? 

We show in Fig. [3] the effect of dropping the 
Ocdm constraint. This is significant in the up- 
per panels, which display the (tooj m i/2) and 
(tan /3, too) planes, but is not so important in 
the (tan/3, TO1/2) and (A), m i/2) planes shown 
in the two lower panels of Fig. [3] These be- 
haviours can be understood by recalling the be- 
haviour of the WMAP coannihilation strips in the 
CMSSM (toi/ 2 ,to ) planes for different values of 
tan (3. For example, the value of too favoured by 
Ocdm for any given values of and Aq in- 

creases as the value of tan j3 increases, foliating 
the {mo,m 1 / 2 ) plane. Thus, for any given value 
of mi/2 and Aq, a large range of values of too can 
be attained for a suitable choice of tan/3, even if 
one does impose the Ocdm constraint. Concern- 
ing the range of m 1 / 2 , this is bounded above by 
(g — 2) and, for any given value of tan/3, the 
allowed range actually decreases for the larger 
values of too allowed if the Ocdm constraint is 
dropped. Thus, dropping the Ocdm constraint 
has little overall effect on the ranges of TO1/2, mo 
and Aq. The primary effect is to enforce a corre- 
lation between tan/3 and toq, as seen in the upper 



right panel of Fig. [3] The range of to decreases 
at any fixed value of tan (3 when the Ocdm con- 
straint is imposed, because of the narrowness of 
the WMAP strip for any fixed value of tan /3. 

Thus, we find that the fit results obtained in 
the parameter planes of too, m\i 2 , Aq and tan/3 
displayed in Fig. [T] are rather robust with respect 
to imposing / dropping the Ocdm constraint. On 
the other hand, as already noted, the Ocdm con- 
straint does reduce the dimensionality of the pa- 
rameter space by essentially one unit. This can 
also be seen from the fact that without impos- 
ing the Ocdm constraint the preferred parameter 
region obtained from the fit to the EWPO and 
BPO still yields a wide range of possible values of 
Ocdm^i 2 - Specifically we find that the 68% C.L. 
region of the (tocTOim) plane shown in Fig. [1] 
yields Ocdm^ 2 < 0.9, while considerably larger 
values of Ocdm^ 2 are allowed at the 95% C.L. 
Eventually, SUSY particle mass measurements at 
the LHC (see the discussion of Fig. [6] below) may 
enable this estimate of Ocdm^ 2 to be refined con- 
siderably (see, e.g., [106|111|112] ). 

Drilling down into the dependences of our 
results on uncertainties in the experimental 
and phenomenological constraints, we display in 
Fig. H] the results of studies of their sensitivi- 
ties to some key observables. The observables 
tracked are (g — 2) M , BR(6 — > S7), Ocdm/i 2 , 
BR(£> U — > tv t ) and mw- The left panel shows 
the percentage variation in the preferred region 
of the (too,to.x/2) plane as the assumed errors in 
these quantities are rescaled, assuming that the 
future experimental central values agree with the 
current ones. The right panel shows the same for 
the area in the (TOo,tan/3) plane. Larger errors 
could arise if we have underestimated the rele- 
vant systematic errors, and smaller errors could 
result from future improvements of the experi- 
mental errors and/or the theoretical predictions. 
As could be expected from the discussion in the 
previous paragraph, the preferred areas vary very 
little with the error in Ocdm, and the areas are 
also relatively insensitive to that in BR(i3„ — > 
tv t ). However, there are greater sensitivities to 
BR(& — > sj), mw and (particularly) (g — 2) M . 

The theoretical error in mw is much smaller 
than the current experimental error. It is en- 
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Figure 3. Variation of the 95% C.L. allowed regions in the MSSM parameter space including (solid) 
or excluding (dotted) the WMAP constraint. The plots show the (mo,m 1 / 2 ) plane (upper left), 
(tan/3, mo) plane (upper right), (tan /3, mx/2) plane (lower left) and the (^4oi m i/2) plane (lower right 
plot). In each panel, we mark the best-fit points found both with and without the WMAP constraint by 
a filled and open star, respectively. 



couraging that reducing the experimental error, 
as should be possible with future Tevatron and 
LHC data, could have substantial effects on the 
preferred areas in the parameter planes. A reduc- 
tion in the error by a factor two could reduce the 



areas by factors of about five, if the present cen- 
tral value (which disagrees with the SM by about 
one a) is maintained. 

The same would be true for a reduction in 
the error in BR(6 — * 57), but here reducing the 
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Figure 4. Relative sizes of the 95% C.L. areas in the (too,TO!/ 2 ) plane (left) and in the (too, tan j3) plane 
(right) as a function of the hypothetical errors of (g — 2) M , BR(6 — > sy), ^cdm/i 2 , BR(_B u — > tv t ), m\y- 
The error scaling is relative to the current combined theory and experimental error. 



theoretical error would also be necessary. This 
would require, in particular, a better understand- 
ing of the uncertainties in higher-order and non- 
perturbative QCD corrections. Indeed, a very 
conservative approach to the combination of the 
current theoretical and experimental errors in 
BR(fo — > sj) might even motivate a larger error 
and hence larger preferred areas than in our de- 
fault analysis. 

Fig. [4] also shows that varying the error in 
(g — 2) M is potentially more important, particu- 
larly if the present error is underestimated. This 
might be the case if, e.g., the weight of experimen- 
tal evidence would shift towards using r decay 
data to estimate the SM hadronic contribution to 
(g — 2)^, or if the error in the light-by-light contri- 
bution were to be revised drastically^. The rapid 
increases in the areas of the preferred regions re- 
flect the fact that a more relaxed treatment of 
the (g — 2) M error led in the past to (parts) of 
the focus-point strip at large too being included 
within the preferred region, which does not occur 
in our default analysis. 

6 In 1 131 it has recently been claimed that solving the muon 
(g — 2) M anomaly by changing the SM prediction of the 
hadronic contribution to (g — 2)^ is unlikely in view of a 
combined analysis of all electroweak data. 



In order to explore the sensitivity to the — 
error in more detail, we show in the left panel of 
Fig. the effect in the CMSSM (to ,to 1/2 ) plane 
of varying this error, while assuming the same 
central value. Going from the outer to the inner 
contours we have assumed erhypothoticai/atoday = 
1.3, 1.2, 1.1, 1.0, 0.9, 0.7 with CTtoday = 8.8 x 10" 10 , 
see Tab. [TJ The partially fuzzy shapes would be 
smoothed by higher statistics. We see that the 
preferred region expands rapidly if the (g — 2) M 
error is increased. Going to an increase by a fac- 
tor of 1.5 (not shown in the plot) would open 
up the focus-point region, which is disfavoured 
in our analysis. Conversely, decreasing this er- 
ror, as would be possible with an accessible im- 
provement of the previous BNL (g — 2) M experi- 
ment [56], would enable the preferred ranges of 
the CMSSM mass parameters to be decreased 
impressively. Ultimately, this together with the 
other EWPO and BPO could make possible a sen- 
sitive test of SUSY at the loop level, if the LHC 
does indeed discover sparticles and measure their 
masses. 

In the right panel of Fig. [5j we make 
a similar analysis of the sensitivity to the 
BR(6 — > S7) error. Going from the outer 
to the inner contours we have again assumed 
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Figure 5. The 95% C.L. region in the (mo> ^1/2) plane for various hypothetical values of the la un- 
certainty (combined theory and experimental) of (g — 2)„(left) and BR(& — > S7)(right). The curves 
show (from outer to inner) the 95% C.L. regions for CThypothcticai/o'today = 1.3,1.2,1.1,1.0,0.9,0.7 for 
fftoday = 8.8 x 10- 10 (left) and CT today (BR^ /BR^ S7 ) = 0.12 (right), respectively see Tab. OQ 



Chypothctical/Otoday — 1.3,1.2,1.1,1.0,0.9,0.7 

with a today (BR^ S7 /BR^ S7 ) = 0.12, sec Tab.ffl 
We see from the right panel of Fig. [5] that treat- 
ing the errors differently could have a significant 
effect. Employing a larger error (as done in [39] . 
for instance), would not only expand the allowed 
regions, but also allow larger tan/3 values, as 
BR(6 — > s-f) is particularly sensitive to tan/3. 

Though in the present analysis we have focused 
on the \i > solution, as favoured by the (g — 2) M 
anomaly, we comment briefly here on the struc- 
ture of the /i < parameter space. In order 
to minimize the discrepancy with the {g — 2)^ 
constraint, for /i < one would need a rela- 
tively heavy spectrum in order to suppress the 
SUSY effects with the wrong sign. This would 
be particularly true for increasing values of tan /3, 
since (g — 2)^ grows almost linearly with tan /3. 
BR(b — > s-f) is also highly sensitive to the sign 
of the /i parameter. In particular, within the 
CMSSM the solution with fi < unambiguously 
implies that all the dominant SUSY effects to 
BR(t> — > s-f) have the same sign and interfere con- 
structively with the SM amplitude. This implies 
more severe constrains with respect to the /i > 



case, and again points toward a heavy spectrum. 
This is not the case for n > 0, where partial can- 
cellations among SUSY effects in BR(fe — > s-f) 
allow relatively light squarks. 

We have also considered possible improvements 
in the determination of the CMSSM parameters 
that might be obtainable from early LHC mea- 
surements. Missing Et measurements with or 
without single lcptons are unlikely to constrain 
the model with high precision. On the other 
hand, in the parameter region preferred by the 
fit (with tan/3 ss 10) there are good prospects 
for measuring the opposite-sign dilepton edge in 
X2 — ► Xi^ + £~ if- = e :^) decays with high preci- 
sion, which is located at 

{m 2 u f^ = (2) 

I 

Such a measurement would constrain a combina- 
tion of sparticle masses and hence the CMSSM 
parameter space in an interesting way. As an ap- 
petizer for what might be possible, we show in 
Fig. [5] the possible impact of a measurement of 
the dilepton edge for the CMSSM best-fit point 
described in the previous paragraph, which has 
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m^o = 121 GeV, m^o = 225 GeV, m r = 
139 GeV, yielding an edge at m e + g - = 87 GeV. 
We assume experimental and theoretical errors of 
3 GeV each. We see in Fig. [5] that the dilepton 
edge measurement would reduce the parameter 
space preferred at the 68% C.L. to two narrow 
strips in the (mo, m 1( / 2 ) plane, linked into a tilted 
'vee' shape at the 95% C.L. The best-fit point in 
the right wing of the 'vee' has quite different pa- 
rameter values from the overall best-fit point in 
the left wing of the 'vee': m 1 / 2 = 390 GeV, mo = 
230 GeV, A = 1230 GeV, tan/3 = 23, yielding 
X 2 = 22.7 and m^o = 155 GeV, m^o = 293 GeV, 
mj = 273 GeV. 

3.3. Comparison with the NUHM1 case 

The above analysis of the CMSSM is relatively 
encouraging for the early days of the LHC, but 
one might wonder to what extent the conclu- 
sions can be extended to more general incarna- 
tions of the MSSM. The full parameter space of 
the MSSM has so many dimensions that exploring 
it with the MCMC approach used here would re- 
quire prohibitive amounts of CPU time. Accord- 
ingly, we discuss briefly here only the simplest 
possible generalization of the CMSSM, in which 
the soft supersymmetry-breaking scalar contribu- 
tions to the Higgs masses are allowed to differ 
by the same amount from those of the squarks 
and sleptons at the GUT scale, the so-called non- 
universal Higgs model 1 (NUHM1) [11311141115] . 

Overall, it is encouraging that the general sizes 
of the 68% and 95% C.L. regions are similar to 
those in the CMSSM, as shown in Fig. though 
the 68% C.L. region together with the best-fit 
point are shifted to lower m.1/2, and the 95% C.L. 
region is more elongated in m^. As in the case of 
the CMSSM, SUSY could be discovered over all of 
the 68% C.L. region with 100 pb _1 of integrated 
luminosity at 14 TeV in a single experiment, and 
even 50 pb _1 of integrated luminosity at 10 TeV 
would cover most of it. As in the CMSSM, not 
all of the NUHM1 95% C.L. region would be cov- 
ered by the LHC with 1 fb _1 of integrated lumi- 
nosity at 14 TeV, whereas the same-sign dilep- 
ton search would cover all the 68% C.L. region in 
the NUHM1. There are differences between the 
shapes of the preferred regions in the CMSSM 



and the NUHM1, particularly at low m.1/2. This 
reflects the fact that the Ocdm constraint can 
be obeyed away from the coannihilation strip at 
larger values of mo, if m x ~ m H /^/2. This free- 
dom can then be exploited to relax the slight ten- 
sion induced by mh which arises in the CMSSM. 

The spectrum at the best-fit NUHM1 point is 
shown in the right panel of Fig. [2] This point 
has m = 100 GeV, m 1/2 = 240 GeV, A Q = 
-930 GeV, tan/3 = 7, m 2 H = -6.9 x 10 5 GeV 2 
and n = 870 GeV, yielding X 2 = 18.0 (39% prob- 
ability) and m/j = 118 GeV. The best- fit values 
of mi/2 and tan (3 are somewhat lower than those 
in the CMSSM, whereas the value of mo is some- 
what higher. The overall value of \ 2 is also some- 
what lower than in the CMSSM, reflecting the re- 
laxation of the slight tension in the value of mh 
that is possible when the Higgs masses are al- 
lowed to become non-universal. Comparing with 
the best-fit CMSSM spectrum, we see that the 
masses of the sleptons and squarks are quite sim- 
ilar, as are the masses of the lighter neutralinos 
and chargino. However, the splitting between the 
stop mass eigenstates is larger — reflecting the 
larger value of \Aq\, the heavier neutralinos and 
chargino are much heavier — reflecting the larger 
value of /x, the best-fit value of mh lies comfort- 
ably above the LEP lower limit, and the heavier 
Higgs bosons are lighter than in the CMSSM — 
reflecting the extra freedom conferred by the non- 
universality in the NUHM1. The lower values of 
the heavier Higgs masses compensate other SUSY 
contributions to BR(6 — » sj), and offer better 
prospects for detection at the LHC than those 
offered by the CMSSM. 

4. Conclusion and Outlook 

Making a probabilistic analysis using a MCMC 
technique, we have presented in this paper the re- 
gions preferred in the CMSSM and the NUHM1 
parameter spaces at the 68% and 95% C.L., as 
well as the spectra at the best-fit points, in the 
light of the present direct and indirect constraints 
on the models' parameters. Particularly impor- 
tant roles are played by (g — 2) M and BK(b — * sj), 
and we have analyzed the ways in which effects of 
these constraints vary with the sizes of their theo- 
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Figure 6. The (too,TOi/2) plane in the CMSSM for tan/3 = 10 and Aq = 0, showing the improvement 
in the constraints on Too and mi/2 that could be obtained by measuring the opposite-sign dilepton edge 
with 1 fb _1 of integrated luminosity at 14 TeV, assuming the spectrum of the best-fit point shown in the 
left panel of Fig. [21 and experimental and theoretical errors of 3 GeV each. The best-fit point is indicated 
by a filled circle. 



retical uncertainties and experimental errors. We 
have quantified how strengthening (or relaxing) 
either of these constraints would reduce (or ex- 
pand) considerably the preferred regions in the 
CMSSM (too, TO^) plane. We have also stud- 
ied the impact of the constraint on the cold dark 
matter density imposed by WMAP. We find that 
the results for the best-fit points are remark- 
ably robust with respect to imposing or drop- 
ping the constraint on the cold dark matter den- 
sity. Perhaps surprisingly, we find that this con- 
straint does not restrict significantly most two- 
dimensional projections of the preferred region in 
the CMSSM parameter space. Encouragingly, we 
find that the preferred regions in the NUHM1 are 
quite similar to those in the CMSSM. 

The 95% exclusion regions in the (too,TOi/2) 
plane extend significantly further than the dis- 
covery regions shown above. Therefore, if SUSY 
were to be excluded at the LHC with 1 fb _1 
(100 pb _1 ) of integrated luminosity at 14 TeV, 
the 95% (68%) C.L. regions in both the CMSSM 



and the NUHM1 would be ruled out. On the 
other hand, SUSY could be discovered at the 5 a 
level at the LHC with 1 fb" 1 of integrated lu- 
minosity at 14 TeV in a single experiment over 
most of the 95% C.L. regions in the (to^to^) 
planes of the CMSSM and the NUHM1. Only 
the highest TO1/2 values would require a larger 
integrated luminosity, or the combination of data 
from both ATLAS and CMS. Indeed, SUSY could 
be discovered over all of the 68% C.L. regions 
in both the CMSSM and the NUHM1 with just 
100 pb -1 of integrated luminosity at 14 TeV, and 
even 50 pb _1 of (good-quality) data at 10 TeV 
would offer significant prospects for SUSY detec- 
tion. The same-sign dilepton search would cover 
most (all) of the 68% C.L. region in the CMSSM 
(NUHM1). If Nature were to choose the best- fit 
CMSSM point, a measurement of the same-sign 
dilepton endpoint would impose a strong con- 
straint on the SUSY spectrum. 

One way or the other, there are good prospects 
that the initial runs of the LHC will determine the 



15 




tan(3= 10, A =0, (i>0 

jets + MET (CMS) 

— ■ leplon + 4 jets (ATLAS) 
• • • 1 lepton + 4 jets (ATLAS) 

- - SS2|j (CMS) 
Higgs (2/fb ) (CMS) 



full CMSSM 
parameter space 

^ 68%C.L. 
t^l 95%C.L. 

NO EWSB 



1000 1200 1400 1600 1800 2000 

m [GeV] 




tan(3= 10, A = 0, \l> 

1/fb@14TeV 

— 100/pb@14TeV 
50/pb@10TeV 





full CMSSM 




parameter space 




^ 68%C.L. 




^ 95%C.L. 



1000 1200 1400 1600 1800 2000 

m [GeV] 



Figure 7. The ("10,^1/2) plane in the CMSSM for tan/3 = 10 and Aq = 0, as in Fig. [TJ overlaid with the 
68% and 95% probability contours for the NUHM1. Upper plot: Some 5 c discovery contours at ATLAS 
and CMS with 1 fb _1 at 14 TeV, and the contour for the 5 c discovery of the Higgs boson in sparticle 
decays with 2 fb _1 at 14 TeV in CMS. Lower plot: The 5cr discovery contours for jet + missing Ej 1 
events at CMS with 1 fb _1 at 14 TeV, 100 pb^ 1 at 14 TeV and 50 pb" 1 at 10 TeV centre-of-mass energy. 



fate of many speculations about the relevance of Acknowledgements 
low-energy SUSY to particle physics. 

We thank A.M. Weber for collaboration in 
the early stages of this work. This work was 
supported in part by the European Commu- 



1G 



nity's Marie-Curie Research Training Network 
under contracts MRTN-CT-2006-035505 'Tools 
and Precision Calculations for Physics Discov- 
eries at Colliders' and MRTN-CT-2006-035482 
'FLAVIAnet', and by the Spanish MEC and 
FEDER under grant FPA2005-01678. The work 
of S.H. was supported in part by CICYT (grant 
FPA 2007-66387), and the work of K.A.O. was 
supported in part by DOE grant DE-FG02- 
94ER-40823 at the University of Minnesota. 



REFERENCES 

1. H. P. Miles, Phys. Rept. 110 (1984) 1. 

2. H. E. Haber and G. L. Kane, Phys. Rept. 117 
(1985) 75. 

3. R. Barbieri, S. Ferrara and C. A. Savoy, Phys. 
Lett. B 119 (1982) 343. 

4. C. Amsler et al. [Particle Data Group], Phys. 
Lett. B 667 (2008) 1. 

5. R. Barate et al. [ALEPH, DELPHI, L3, 
OPAL Collaborations and LEP Working- 
Group for Higgs boson searches], Phys. Lett. 
B 565 (2003) 61 [arXiv:hep-ex/0306033| . 

6. S. Schael et al. [ALEPH, DELPHI, L3, OPAL 
Collaborations and LEP Working Group for 
Higgs boson searches], Eur. Phys. J. C 47 

(2006) 547 [arXiv:hep-ex/0602042| . 

7. G. W. Bennett et al. [Muon G-2 Collab- 
oration], Phys. Rev. D 73 (2006) 072003 
|arXiv:hep-ex/0602035l . 

8. T. Moroi, Phys. Rev. D 53 (1996) 
6565 [Erratum-ibid. D 56 (1997) 4424] 
[arXiv:hep-ph/9512396] . 

9. A. Czarnecki and W. J. Marciano, Phys. Rev. 
D 64 (2001) 013014 [arXiv:hep-ph/0102122l . 

10. M. Davier, Nucl. Phys. Proc. S uppl. 169 

(2007) 288 [arXiv:hep-ph/0701163| . 

11. J. P. Miller, E. de Rafael and B. L. Roberts, 
Rept. Prog. Phys. 70, 795 (2007) 
|arXiv:hep-ph/0703049] . 

12. F. Jege rlehner, Acta Phys. Pol on. B 38, 3021 
(2007) [arXiv:hep-ph/0703125| . 

13. M. Passera, W. J. Marciano and A. Sir- 
lin, Phys. Rev . D 78, 013009 (2008) 
|arXiv:0804.lT42l [hep-ph]]. 

14. J. Dunkley et al. [WMAP Collaboration], 
larXiv:0803.0586l [astro-ph], 

15. B. C. Allanach and C. G. Lester, Corn- 
put. Phys. Commun. 179 (2008) 256 
|arXiv:0705.0486l [hep-ph]]. 

16. A. Djouadi, M. D rees and J. L. Kneur, J HEP 
0108 (2001) 055 [arXiv:hep-ph/0107316| . 

17. W. de Boer, M. Huber, C. Sander and 
D. I. Kazakov, Phys. Lett. B 515 (2001) 283. 

18. W. de Boer and C. Sander, Phys. L ett. B 585 
(2004) 276 [arXiv:hep-ph/0307049| . 

19. J. R. Ellis, K. A. Olive and Y. San- 
toso, New Jour. Phys. 4 (2002) 32 



17 



[arXiv:hep-ph/0202110| . 

20. H. Baer, C. Balazs, A. Belyaev, 
J. K. Mizukoshi, X. Tata and Y. Wang, JHEP 
0207 (2002) 050 |arXiv:hep-ph/0205325| . 

21. R. Arnowitt and B. Dutta, 
arXiv:hep-ph/0211417 

22. J. R. Ellis, K. A. Olive, Y. Santoso and 
V. C. Spanos, Phys. Le tt. B 565 (2003) 176 
[arXiv:hep-ph/0303043] . 

23. H. Bae r and C. Balazs, JCA P 0305, 006 

(2003) [arXiv:hep-ph/0303114] . 

24. A. B. Lahanas and D. V. Nanopou- 
los, Phys. Lett. B 568, 55 (2003) 
|arXiv:hep-ph/0303130] . 

25. U. Chattopadhyay, A. Corsetti and 
P. Nath, Phys. Rev. D 68, 035005 (2003) 
|arXiv:hep-ph/030320l] . 

26. C. Mu noz, Int. J. Mod. Phy s. A 19, 3093 

(2004) [arXiv:hep-ph/0309346l . 

27. G. Belanger, F. Boudjema, A. Cottrant, 
A. Pukhov and A. Semenov, Nucl. Phys. B 
706 (2005) 411 [arXiv:hep-ph/0407218| . 

28. J. R. Ellis, K. A. Olive, Y. Santoso and 
V. C. Spanos, Phys. Rev . D 69 (2004) 095004 
|arXiv:hep-ph/0310356] . 

29. J. R. Ellis, S. Heinemeyer, K. A. Olive 
and G. Weiglein, JHE P 0502 (2005) 013 
|arXiv:hep-ph/0411216] . 

30. J. R. Ellis, D. V. Nanopoulos, K. A. Olive 
and Y. Santoso, Phys. L ett. B 633 (2006) 583 
|arXiv:hep-ph/050933l] . 

31. J. R. Ellis, S. Heinemeyer, K. A. Olive 
and G. Weiglein, JHE P 0605 (2006) 005 
[arXiv:hep-ph/0602220] . 

32. J. Ellis, S. Heinemeyer, K.A. Olive, A.M. We- 
ber, G. Weiglei n, JHEP 08 (2007) 083 
[arXiv:0706.0652l [hep-ph]] 

33. J. R. Ellis, S. Heinemeyer, K. A. Olive and 
G. Weiglein, Phys. Lett. B 653 (2007) 292 
[arXiv:0706.0977l [hep-ph]]. 

34. E. A. Balt z and P. Gondolo, JH EP 0410 
(2004) 052 [arXiv:hep-ph/0407039] . 

35. B. C. Allanach and C. G. Lester, Phys. Rev. 
D 73 (2006) 015013 [arXiv:hep-ph/0507283| . 

36. B. C. Allanach, Phys. L ett. B 635 (2006) 123 
arXiv:hep-ph/0601089l . 

37. B. C. Allanach, C. G. Lester and 
A. M. Weber, JHEP 0612 (2006) 065 



|arXiv:hep-ph/0609295l . 

38. B. C. Allanach, C. G . Lester and A. M. We- 
ber, larXiv:0705.0487l [hep-ph] . 

39. F. Feroz, B. C. Allanach, M. Hobson, 
S. S. AbdusSalam, R. Trotta and A. M. We- 
ber, larXiv:0807.45l2l [hep-ph] . 

40. R. R. de Austri, R. Trotta and 
L. Roszkowski, JHEP 0605 (2006) 002 
|arXiv:hep-ph/0602028| . 

41. L. Roszkowski, R. R. de Austri and 
R. Trotta, JHEP 0704 (2007) 084 
[arXiv:hep-ph/0611173| . 

42. L. Roszkowski, R. Ruiz de Austri and 
R. Trotta, JH EP 0707 (2007) 075 
[arXiv:0705.20T2l [hep-ph]]. 

43. L. Roszkowski, R. R. de Austri, J. Silk and 
R. Trotta. larX~iv:0707.0622l [astro-ph], 

44. S. Heinemeyer, X. Miao, S. Su and 
G. Weiglein, JHEP 0808 (2008) 087 
[arXiv:0805.2359l [hep-ph]]. 

45. P. Bechtle, K. Desch and P. Wienemann, 
Comput. Phys. Comm un. 174 (2006) 47 
[arXiv:hep-ph/04120T2| . 

46. R. Lafaye, T. Plehn, M. Rauch and D. Zer- 
was, Eur. Phy s. J. C 54 (2008) 617 
[arXiv:0709.3985l [hep-ph]]. 

47. O. Buchmueller et al, Phys. Lett. B 657 
(2007) 87 |arXiv:0707.3447l [hep-ph]]. 

48. B. C. Allanach and D. Hooper, 
larXiv:0806.1923l [hep-ph] . 

49. R. Trotta, Talk given at tools08, 
http : / / indi co"! cern. ch/getFile -py/| 
access?contribld=43&sessionld=20& 
resld=0&materialld=slides&conf Id=35476, 
based on work done in collaboration with 
L. Roszkowski and R. R. de Austri. 

50. S. Heinemeyer, W. Hollik, D. Stockinger, 
A. M. Web er and G. Weiglein, JHEP 0608 
(2006) 052 [arXiv:hep-ph/0604147] . 

51. S. Heinemeyer, W. Hollik, A. M. Weber 
and G. Weiglein , JHEP 0804 (2008) 039 
[arXiv:0710.2972l [hep-ph]] . 

52. M. Verzocchi, talk at ICHEP 2008, August 
2008, Philadelphia, USA. 

53. G. Degrassi and G. F. Giudice, Phys. Rev. D 
58 (1998) 053007 [arXiv:hep-ph/9803384l . 

54. S. Heinemeyer, D. Stockinger and G. Wei- 
glein, Nucl. Phys. B 690 (2004) 62 



18 



|arXiv:hep-ph/0312264| . 

55. S. Heinemeyer, D. Stockinger and G. Wei- 
glein, Nucl. Phys. B 699 (2004) 103 
|arXiv:hep-ph/0405255l . 

56. D. W. Hertzog, J. P. Miller, E. de 
Rafael, B. Lee Roberts and D. Stockinger, 
larXiv:0705.46i~7l [hep-ph] . 

57. G. Degrassi, S. Heinemeyer, W. Hollik, 
P. Slavich and G. Weiglein, Eur. Phy s. J. C 
28 (2003) 133 [arXiv:hep-ph/0212020| . 

58. S. Heinemeyer, W. Hollik and G. Wei- 
glein, Eur. Phys. J . C 9 (1999) 343 
|arXiv:hep-ph/9812472| . 

59. S. Heinemeyer, W. Hollik and G. Wei- 
glein, Comput. Phys. Commun. 124 
(2000) 76 [arXiv:hep-ph/9812320| . See 
|http : //www . f eynhiggs . de 

60. M. Frank, T. Hahn, S. Heinemeyer, W. Hol- 
lik, H. Rzeh ak and G. Weiglein, JH EP 0702 
(2007) 047 [arXiv:hep-ph/0611326| . 

61. M. Misia k et al., Phys. Rev. L ett. 98 (2007) 
022002 |arXiv:hep-ph/0609232| . 

62. M. Ciuchini, G. Degrassi, P. Gambino and 
G. F. Giudice, Nucl. P hys. B 534 (1998) 3 
arXiv:hep-ph/9806308l . 

63. G. Degrassi, P. Gambino and 
G. F. Giudice, JHEP 0012 (2000) 009 
arXiv:hep-ph/0009337l . 

64. M. S. Carena, D. Garcia, U. Nierste and 
C. E M. Wagner, Phys. Le tt. B 499 (2001) 
141 |arXiv:hep-ph/0010003| . 

65. G. D'Ambrosio, G. F. Giudice, G. Isidori and 
A. Strumia, Nucl. Phy s. B 645 (2002) 155 
|arXiv:hep-ph/0207036| . 

66. E. Barberio et al. [Heavy Flavour Av- 
eraging Group (HFAG)] , |hep-ex/0603003| 
http : / / slac . Stanford . edu/xorg/hf ag/ . 

67. Tevatron Electroweak Working Group, for the 
CDF and DO Collaborations, [MXiyi)80jUJ)89J 
[hep-ex] . 

68. G. Belanger, F. Boudjema, A. Pukhov and 
A. Semenov, Comput. Phys. Commun. 176 
(2007) 367 [arXiv:hep-ph/0607059| . 

69. G. Belanger, F. Boudjema, A. Pukhov and 
A. Semenov, Comput. Phys. Commun. 149 
(2002) 103 [arXiv:hep-ph/0112278| . 

70. G. Belanger, F. Boudjema, A. Pukhov and 
A. Semenov, Comput. Phys. Commun. 174 



(2006) 577 [arXiv:hep-ph7 0405253|. 

71. G. Isidori and A. Retico, JH EP 0111 (2001) 
001 |arXiv:hep-ph/0110121| . 

72. A. J. Buras, P. H. Chankowski, J. Rosiek and 
L. Slawianowska, Nucl. Phys. B 659 (2003) 3 
[hep-ph/0210145] ; 

73. G. Isidori a nd P. Paradisi, Phys. L ett. B 639 
(2006) 499 [arXiv:hep-ph/0605012| . 

74. G. Isidori, F. Mescia, P. Paradisi and 
D. Te mes, Phys. Rev. D 75 (2007) 
115019 |arXiv:hep-ph/0703035] , and refer- 
ences therein. 

75. A. G. Akeroyd and S. Recksieg el, J. Ph ys. G 
29 (2003) 2311 [arXiv:hep-ph/0306037l . 

76. B. Aubert et al. [BABAR Collabora- 
tion], Phys. Rev. L ett. 95 (2005) 041804 
|arXiv:hep-ex/0407038| . 

77. P. Chang, talk at ICHEP 2008, August 2008, 
Philadelphia, USA. 

78. A. Gray et al. [HPQCD Collaboration], 
Phys. Rev. L ett. 95 (2005) 212001 
[hep-lat/0507015| . 

79. C. Bobeth, A. J. Buras and T. Ew- 
erth, Nucl. Phys. B 713 (2005) 522 
|arXiv:hep-ph/0409293] . 

80. T. Huber, E. Lunghi, M. Misiak and 
D. Wyler, Nucl. Ph ys. B 740 (2006) 105 
|arXiv:hep-ph/0512066] . 

81. M. Antonelli et al. [FlaviaNet Working Group 
on Kaon Decays]. larXiv:0801. 18171 [hep-ph], 

82. A. J. Buras, P. Gambino, M. Gorbahn, 
S. Jager a nd L. Silvestrini, Nucl. P hys. B 592 
(2001) 55 [arXiv:hep-ph/0007313| . 

83. A. V. Artamonov et al. [The E949 Collabora- 
tion] , larXiv:0808.2459l [hep-ex] . 

84. M. Bona et al. [UTfit Collaboration], JHEP 
0803 (2008) 049 [arXiv:0707.0636l [hep-ph]] . 

85. T. Becher and M. Neubert, Phys. Rev. L ett. 
98 (2007) 022003 [arXiv:hep-ph/0610067| . 

86. P. Gambino and P. Giordano, 
larXiv:0805.027ll [hep-ph]. 

87. Z. Ligeti, I. W. Stewart and F. J. Tackmann, 
larXiv:0807.l926l [hep-ph] . 

88. M. Misiak. [arXiv:0808.3134l [hep-ph], 

89. P. Gambino, U. Haisch and M. Misiak, 
Phys. Rev. Lett. 9 4, 061803 (2005) 
|arXiv:hep-ph/0410155] . 

90. T. Hurth, G. Isidori, J. F. Kamenik, F. Mes- 



19 



cia, iarXiv:0807.50 39 [hep-ph]. 

91. C. Bobeth, M. Misiak and J. Ur- 
ban, Nucl. Phys. B 567, 153 (2000) 
[arXiv:hep-ph/9904413] . 

92. G. Degrassi, P. Gambino and P. Slavich, 
larXiv:0712.3265l [hep-ph] . 

93. F. Mahmoudi, Comput. Phys. Commun. 178 
(2008) 745 [arXiv:0710^2067l [hep-ph]] and 
larXiv:0808.31~44l [hep-ph]. 

94. J. R. Ellis, S. Heinemeyer, K. A. Olive and 
G. Weiglein, Phys. Let t. B 515 (2001) 348 
arXiv:hep-ph/010506l] . 

95. S. Ambrosanio, A. Dedes, S. Heinemeyer, 
S. Su an d G. Weiglein, Nucl. P hys. B 624 
(2002) 3 [arXiv:hep-ph/0106255| . 

96. B. C. Allanach, Comput. Phys. Commun. 143 
(2002) 305 |arXiv:hep ph/0104145| . 

97. P. Gondolo, J. Edsjo, P. Ullio, L. Bergstrom, 
M. Schelke and E. A. Baltz, New Astron. Rev. 
49 (2005) 149. 

98. P. Gondolo, J. Edsjo, P. Ullio, L. Bergstrom, 
M. Schelke and E. A. Baltz, JCA P 0407 
(2004) 008 [arXiv:astro-ph/0406204] . 

99. P. Skands et al, JH EP 0407 (2004) 036 
[arXiv:hep-ph/031 1 123] . 

100. F. James and M. Roos, Comput. Phys. Com- 
mun. 10 (1975) 343. 

101. ATLAS Collaboration, Detector and Physics 
Performance Technical Design Report, 
CERN/LH CC/99-15 (1999), see: 

http : //atlasinf o . cern. ch/Atlas/ 
GROUPS/PHYSICS/TDR/access.html ; 
https : / / twiki . cern . ch/ twiki/bin/ view/ 
Atlas/SUSYWorkingGroup . 

102. G. L. Bayatian et al, CMS Collab- 
oration, CMS Technical Design Re- 
port, Volume II: Physics Performance, 
CERN-LHCC-2006-021, CMS-TDR- 
008-2 J. Phys. G34, 995 (2007); see: 
|http : / / cmsdoc . cern . ch/ cms/ cpt/tdr/ . 

103. J.-J. Blaising, A. De Roeck, J. Ellis, F. Gian- 
otti, P. Janot, G. Rolandi, and D. Schlatter, 
Potential LHC Contributions to Europe 's Fu- 
ture Strategy at the High Energy Frontier, 

J. Ellis, |arXiv:hep-ph/0611237| 

104. B. A llanach et al, Eur. Phy s. J. C 25 (2002) 
113 [arXiv:hep ph/0202233] . 

105. M. Battaglia et al, Eur. Phys. J. C 22 (2001) 



535 [arXiv:hep-ph/0106204] . 

106. M. Battaglia, A. De Roeck, J. R. Ellis, F. Gi- 
anotti, K. A. Olive and L. Pape, Eur. Phys. 
J. C 33 (2004) 273 |arXiv:hep-ph7 0306219 . 

107. G. Weiglein et al [LHC/LC Study 
Group], Phys. Re pt. 426 (2006) 47 
[arXiv:hep-ph/0410364| . 

108. T. Plehn, |http : //www . ph . ed . ac . uk/ 
^tplehn/prospino/ . 

109. V. M. Abazov et al [DO Collaboration], Phys. 
Lett. B 660 (2008) 449 [arXiv:0712. 3805 [hep- 
ex]]. 

110. S. Dube et al. [CDF Collaboration], CDF 
public note 9176, 

http : //www-cdf . f nal . gov/physics/ 
exotic/r2a/20080110.trileptonjdube . 

111. M. M. Nojiri, G. Polesello and 
D. R. Tovey, JHEP 0603 (2006) 063 
|arXiv:hep ph/0512204] . 

112. E. A. Baltz, M. Battaglia, M. E. Peskin and 
T. Wizansky, Phys. Rev . D 74 (2006) 103521 
|arXiv:hep-ph/0602 187| . 

113. H. Baer, A. Mustafayev, S. Profumo, 
A. Belyaev an d X. Tata, Phys. Rev. D 71, 
095008 (2005) |arXiv:hep-ph/0412059| . 

114. H. Baer, A. Mustafayev, S. Profumo, 
A. Belyaev and X. Tata, J HEP 0507 (2005) 
065 |arXiv:hep-ph/0504001| . 

115. J. Ellis, K. A. Olive and P. Sandick, 
larXiv:0805.2343l [hep-ph] . 



